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ABSTRACT 

We study the rest-frame ultra-violet sizes of massive (~ 0.8 x 10 ®®Mq) galaxies at 3.4 < z < 4.2, 
selected from the FourStar Galaxy Evolution Survey (ZFOURGE), by fitting single Sersic profiles 
to HST/WFG3/FI60W images from the Gosmic Assembly Near-Infrared Deep Extragalactic Legacy 
Survey (GANDELS). Massive quiescent galaxies are very compact, with a median circularized half- 
light radius re = 0.63 ± O.ISkpc. Removing 5/16 (31%) sources with signs of AGN activity does not 
change the result. Star-forming galaxies have re = 2.0 ± 0.60kpc, 3.2 ± I.3x larger than quiescent 
galaxies. Quiescent galaxies at z 4 are on average 6.0 ± I.7x smaller than a.t z ^ 0 and 1.9 ± 0.7x 
smaller than at z ^ 2. Star-forming galaxies of the same stellar mass are 2.4 ± 0.7x smaller than 
at z 0. Overall, the size evolution at 0 < z < 4 is well described W a powerlaw, with re = 

5.08 ± 0.28(1 -k z)-®'44±° °®kpc for quiescent and = 6.02 ± 0.28(1 -k z)-°'^2±o.05]jpp star-forming 
galaxies. Compact star-forming galaxies are rare in our sample: we find only 1/14 ^ 7% with 
re/(M/lO®®M 0 )°-®’^ < 1.5, whereas 13/16 ^ 81% of the quiescent galaxies is compact. The number 
density of compact quiescent galaxies at z ^ 4 is 1.8 ± 0.8 x 10“®Mpc“^ and increases rapidly, by 
> 5x, between 2 < z < 4. The paucity of compact star-forming galaxies at z ^ 4 and their large 
rest-frame ultra-violet median sizes suggest that the formation phase of compact cores is very short 
and/or highly dust obscured. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — infrared: 
galaxies — cosmology: observations 


1. INTRODUCTION 


In recent years massive qu iescent galaxies have 
been found beyond z = 3 (e.g. iChen fc Marzkel 12004 
Wiklind et al.l 120 081: IMancini et al.l 120091 F ontana et al l 


2009t^_ MarcheshhetaL _ 20101 jGuoe^^iy__l201^ 

Stefanon et al .1 l2013t iMuzzin et al.l 120131 iSnitler et iil] 


201411 and even at z ^ 4, when the universe was 


only 1.5 Gyrs old (|Straatman et al.l l2014f) . Quiescent 
galaxies at high redshift (z > 1 ) exhi bit compact 
morphologies, with small effective radii (e.g[DaddLet_alJ 
120051 : Ivan Dokkiim et aP 120081 : iDamianov et al.l 12009(1 . 
which tend to become smaller with increasing redshift 


^ straatman@strw.leidenuniv.nl 

^ Leiden Observatory, Leiden University, PO Box 9513, 2300 
RA Leiden, The Netherlands 

^ Australian Astronomical Observatory, PO Box 915, North 
Ryde, NSW 1670, Australia 

Department of Physics Astronomy, Macquarie University, 
Sydney, NSW 2109, Australia 

^ Centre for Astrophysics and Supercomputing, Swinburne 
University, Hawthorn, VIC 3122, Australia 

® George P. and Cynthia W. Mitchell Institute for Fundamen¬ 
tal Physics and Astronomy, Department of Physics and Astron- 
orrw, Texas A&; M University, College Station, TX 77843 

' Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD 21218, USA 

® Department of Astronomy, Yale University, New Haven, CT 
06520, USA 

® Carnegie Observatories, Pasadena, CA 91101, USA 

Astrophysics Science Division, Goddard, Space Flight Cen¬ 
ter, Code 665, Greenbelt MD 20771, USA 

This paper contains data gathered with the 6.5 meter Mag¬ 
ellan Telescopes located at Las Campanas observatory, Chile. 


(|van der Wei et al'1l2014[ f. At z 3, they have sizes of 
^ Ikpc, 3 —4x smaller t han early-type galaxies of similar 
stellar mass at z ~ 0 (|Shen et al.l 120031 : iMosleh et al.l 
120131) and 2 — 3x smaller than star-forming galaxies at 
the same redshift. 

How compact quiescent galaxies are formed is still un¬ 
clear. Simulations propose mechanisms in which gas- 
rich major mergers can induce cent ral starbursts. result¬ 
ing in a compact m erger remnant (jHoDkins et al.l [20091 : 

I Wellons et al.l [ 20 Tll . or in which massive star-forming 
clu mps move to the c enters if galaxy disks are unsta¬ 
ble (iDekel et al.l 1200^ iDekel fc BurkerO 120141) . Alterna¬ 
tively they may have formed in a more protracted pro- 
cess at high red shift. when the universe was more dense 
(|Mo et al.lll998D . 

To understand these scenarios it is necessary to identify 
compact quiescent galaxies and their progenitors at the 
highest redshifts. Compact star-f orming galaxies be en 
found in sma l l num bers at z = 2—3 (jBarro et al.ll2014al lbl: 
iNelson et al.ll20fl) . but many host AGN, complicating 
the interpretation of the observations. At the same time, 
rest-frame ultra-violet (UV) or optically measured sizes 
of star-forming galaxies may be affected by dust-obscured 
central regions, thereby increasing their effective radii. 

In this work we investigate the sizes of a stellar-mass 
complete sample of star-forming and quiescent galax¬ 
ies at z ^ 4. Throughout, we assume a standard 
ACDM cosmology with Dm = 0.3, Da = 0.7 and Hq = 
70km s“®Mpc“®. The adopted photometric system is 
AB. 
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2. SAMPLE SELECTION 

The galaxies were selected using deep iA^-band 
images from the FourStar Galaxy Evolution Survey 
(ZFOURGE; Labbe et al. in prep .), a near-IR surve y 
with the FourStar Infrared Camera (jPersson et al.l[2M^ . 
coverin g three 11 ' x 11 ' poin tings, locate d in the fields 
CPF S (jGiaccon i et al.||2 002D. COSM OS (jScoville et al.l 
[2?in7ll and UPS dbawrence et al1[2?i?il . The ZFOURGE 
A"s-band selected catalogs are at least 80% complete 
down to Kg = 24.53 , 24.74 and 25.07 mag in 
each field, respectively (iPapovich et al.l I2014D . Pho¬ 
tometric redshifts and stellar masses were derived us¬ 
ing 5 near-IR medium-bandwidth filters on FourStar 
{Ji, J 2 , J 3 , Hg, Hi), which provide fine sampling of the 
age-sensitive Balmer/4000A break at 1.5 < z < 

4, in combination wit h public data over a wave¬ 
length range 0.3 — 8 /nn (iStraatman et al.ll2014D . Here 
we make additio nal use of HST/WFC3/F160W data 
from CANPELS (iGrogin et all 120111 iKoekemoer et all 
120111 iSkelton et al.l 1201411 . to examine galaxy sizes and 
Spitzer/MIPS 24^m data from GOOPS-South (PI: Pick- 
inson), COSMOS (PI: Scoville) and SPUPS (PI: Punlop) 
to measure infrared flux. 

The galaxies in this work have photometric redshifts 
3.4 < 2 : < 4.2, stellar masses of logiQ(M/M 0 ) > 
10.55 and ATg-band signal-to-noise (SNR) of SNR> 
7. They are separated into quiescent and star¬ 
forming accord ing to their rest-frame U — V versus 
y — J colours (iLabbe et ^ 120051 iWilliams et al.l 120091 
iSpitler et al.l 120141) yielding 19 quiescent and 25 star¬ 
forming galaxies (jStraatman et al.l I2014H . Of these, 
34 have HST/WFC3/F160W coverage. One quiescent 
galaxy has SNR< 3 in F160W and is not included. An¬ 
other star-forming galaxy with a highly uncertain red- 
shift solution was also rejected from the sample, along 
with two star-forming galaxies that appear to consist of 
two sources each in the higher resolution HST images. In 
total we study 16 quiescent and 14 star-forming galaxies. 
We include a control sample at 2 < 2 < 3.4 (326 sources) 
at similar mass and SNR. 

3. GALAXY SIZES FROM HST/WFC3 IMAGING 
3.1. Sersic fits 

Sizes and struc tural parameters were measured 
by fitting Sersic (|Sersi 3 I1968H profiles on 6 " X 6 " 
HST/WFC3/F16 0W image stamps using GALFIT 
(|Peng et al.ll201(i|l . In particular, we measure the half- 
light radius, encapsulating half the sources’ integrated 
light. The corresponding parameter in GALFIT is the 
half-light radius along the semi-major axis {ri/ 2 ^maj)-i 
which can be converted to circularized effective radius 
(re = ri/ 2 ,maj\/{b/a)), with h/a the axis ratio. 

We manually subtracted the background in each image 
stamp, masking sources and using the mode of the pixel 
flux distribution. Sky estimation in GALFIT was turned 
off. Neighbouring objects at r > 1.1" from the source 
were effectively masked by setting their corresponding 
pixels in the image to zero flux and increasing those in 
the noise image by xlO®. Glose neighbouring objects 
were fitted simultaneously. 

We created mean PSFs for each field by stacking im¬ 
age stamps of bright stars (masking all neighbouring 
sources). As many of the galaxies are small we inves¬ 


tigate the impact of PSF choic e. We repeated the fitting 
using the hybrid PSF models of Ivan der Wei et all (| 201 ^ 
and find marginally larger (< 5%) sizes. In particular, 
for the smallest galaxies (re < 0 . 20 "), we find a median 
re/'f’e,PSFvdW = 0.93 ± 0.05. 

Errors on the individual measurements were calculated 
using a Monte Carlo procedure. After subtracting the 
best-fit GALFIT models from the sources, we shifted the 
residuals by a random number of pixels, added back the 
model and used this as input for GALFIT. Repeating 
this > 200 X for each galaxy, errors were calculated as 
the 1(7 variation on these measurements. We report our 
results in Tabled] 

In the fits, the Sersic index (ngersic) was restricted to 
0.1 < Ugersic < 8.0. If Ugersic reached the extreme value 
0.1 or 8.0, GALFIT was rerun while forcing Ugersic = 1 
for star-forming and ng^rsic = 4 for quiescent galaxies. 
These values correspond to the median Ugersic of galaxies 
with well-constrained fits and SNRfisow >15. 

At 2 ~ 4, this happens for 6/16 => 38% quies¬ 
cent and 2/14 14% star-forming galaxies. To ex¬ 

plore systematic effects introduced by the choice of pro¬ 
file, we set Ugersic = 1-0 or Ugersic = 4.0 for bright 
{magpieowi^^B) < 24.5) and compact sources {re < 
0.20") and find on average re,n=i/?'e,n =4 = 0.80 ± 0.13, 
corresponding to a systematic uncertainty of 20%. We 
add this in quadrature to the uncertainties from the 
Monte Garlo procedure for each galaxy. Systematic bi¬ 
ases of t his level do not affe c t the m ain results. For com¬ 
parison, Ivan der Wei et all (120121) derived typical sys¬ 
tematic uncertainties on size of ~ 12% for faint F160W= 
24 — 26 and small Ve < 0.3" galaxies. 

As many galaxies have modest SNR, we tested the re¬ 
liability of our measurements by a simulation, in which 
we inserted source models, convolved with the instru¬ 
ment PSF, in the F160W images. These have adopted 
magnitudes of 25 < ruagpieowi^B) < 26 and size of 
0.06 < re(") < 0.3. We find re,outlre,in = 0.97 ± 0.05, 
with Te^in and re,out the input and output effective radii, 
showing that we can recover the sizes of faint compact 
sources without bias. As an additional test we deter¬ 
mine the size distribution of point sources, by inserting 
PSFs in the images and measuring their size. We can 
constrain the size of bright objects to 0.01" at 95% con¬ 
fidence, which we adopt as a minimum uncertainty on 
the sizes. 

We crossmatc hed our sample at 2 < 2 < 4.2 with the 
size catalogs of Ivan der Wei et"^ (I2014D. based o n the 
3D—HST photometric catalogs (|Skelton et al.ll2014l) . We 
find that the sizes and Sersic indices agree well, with 
a median re,zFOURGE/re, 3 DHST = 1.004 ± 0.01 and 
UZFOURGE — UsdHST = —0.012 ± 0.058. 

We test for color gradients between rest-frame UV sizes 
and rest-frame optical sizes, using a rest-frame color and 
stellar-mass matched control sample at 2 3. We find 

F160W (rest-frame 4000A) sizes are 0 ± 6 % and 6 ± 11 % 
smaller than F125W (rest-frame 3000A) sizes for star¬ 
forming and quiescent galaxies, respectively. 

3.2. Stacking 

We also measure the average sizes by stacking the back¬ 
ground subtracted image stamps of the two subsamples, 
normalizing each by mean stellar mass. Neighbouring 
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QUIESCENT 


ZP-CDPS-209 

53.1132774 

-27.8698730 

3.56 

22.6 

24.1 

24.3 

± 

0.0 

64.6 

0.76 

0.23 

0.06 

± 

0.01 

0.27 

± 

0.07 

0.37 

± 

0.08 

4.00 



0.3 

-0.9 ± 3.5 

ZP-CDPS-403 

53.0784111 

-27.8598385 

3.660*^ 

22.4 

23.7 

23.5 

± 

0.0 

118.0 

1.15 

0.22 

0.12 

± 

0.03 

0.82 

± 

0.18 

0.85 

± 

0.05 

7.78 

± 

0.94 

0.8 

99.8 ± 148.5t 

ZP-CDPS-4719 

53.1969414 

-27.7604313 

3.59 

23.4 

25.2 

25.2 

± 

0.1 

33.5 

0.45 

0.23 

0.12 

± 

0.03 

0.60 

± 

0.14 

0.48 

± 

0.08 

1.88 

± 

0.84 

0.3 

1.9 ± 3.4 

ZP-CDPS-4907 

53.1812820 

-27.7564163 

3.46 

23.6 

25.0 

25.1 

± 

0.1 

38.2 

0.40 

0.28 

0.08 

± 

0.02 

0.56 

± 

0.13 

0.86 

± 

0.12 

3.28 

± 

0.90 

0.8 

1.4 ± 3.6 

ZP-CDPS-5657 

53.0106506 

-27.7416019 

3.56 

23.0 

24.6 

24.2 

± 

0.1 

26.7 

0.76 

0.33 

0.52 

± 

0.16 

3.22 

± 

0.93 

0.72 

± 

0.11 

4.45 

± 

0.98 

0.3 

1.7 ± 3.8t 

ZP-CDPS-617 

53.1243553 

-27.8516121 

3.700*^ 

22.3 

23.5 

23.5 

± 

0.0 

135.1 

0.69 

0.22 

0.10 

± 

0.02 

0.55 

± 

0.11 

0.59 

± 

0.03 

4.00 



0.3 

86.3 ± 3.4t* 

ZP-COSMOS-13129 

150.1125641 

2.3765368 

3.81 

23.6 

25.2 

24.9 

± 

0.1 

10.8 

1.78 

0.46 

0.52 

± 

0.13 

2.15 

± 

0.48 

0.34 

± 

0.08 

0.56 

± 

0.24 

0.6 

110.1 ± 10.2* 

ZP-COSMOS-13172 

150.0615082 

2.3786869 

3.55 

22.4 

24.4 

24.4 

± 

0.1 

37.2 

1.45 

0.27 

0.08 

± 

0.02 

0.49 

± 

0.12 

0.64 

± 

0.13 

3.94 

± 

1.11 

0.6 

2.7 ib 7.6 

ZP-COSMOS-13414 

150.0667114 

2.3823516 

3.57 

23.4 

25.4 

25.4 

± 

0.1 

14.0 

0.44 

0.32 

0.20 

± 

0.06 

0.83 

± 

0.29 

0.34 

± 

0.14 

1.51 

± 

1.00 

0.2 

7.1 ± 8.7 

ZP-UDS-10684 

34.3650742 

-5.1488328 

3.95 

24.1 

25.9 

25.2 

± 

0.2 

8.5 

0.85 

0.32 

0.50 

± 

0.17 

2.42 

± 

0.77 

0.47 

± 

0.18 

4.63 

± 

1.68 

1.0 

8.8 ± 12.8 

ZP-UDS-11483 

34.3996315 

-5.1363320 

3.63 

23.6 

26.0 

25.9 

± 

0.2 

8.9 

1.02 

0.35 

0.11 

± 

0.05 

0.52 

± 

0.25 

0.43 

± 

0.24 

4.59 

± 

2.01 

1.0 

1.8 ± 10.2 

ZP-UDS-2622 

34.2894516 

-5.2698011 

3.77 

23.0 

24.6 

24.5 

± 

0.1 

29.9 

0.87 

0.30 

0.13 

± 

0.03 

0.76 

± 

0.19 

0.66 

± 

0.10 

4.00 



0.9 

12.2 ± 10.6 

ZP-UDS-3112 

34.2904282 

-5.2620673 

3.53 

23.2 

24.9 

24.9 

± 

0.1 

25.7 

0.43 

0.30 

0.07 

± 

0.02 

0.39 

± 

0.13 

0.66 

± 

0.19 

4.00 



0.0 

-10.9 ± 10.6 

ZP-UDS-5418 

34.2937546 

-5.2269468 

3.53 

23.3 

24.9 

24.9 

± 

0.1 

20.7 

0.44 

0.30 

0.07 

± 

0.02 

0.50 

± 

0.14 

0.83 

± 

0.17 

4.00 



0.5 

48.4 ± 10.6 

ZP-UDS-6119 

34.2805405 

-5.2171388 

4.05 

23.8 

25.5 

25.4 

± 

0.2 

10.6 

0.55 

0.32 

0.26 

± 

0.15 

1.26 

± 

0.75 

0.49 

± 

0.20 

4.00 



1.0 

-12.5 ± 8.7 

ZP-UDS-9526 

34.3381844 

-5.1661916 

3.97 

24.2 

25.9 

25.8 

± 

0.3 

11.5 

0.89 

0.21 

0.10 

± 

0.05 

0.39 

± 

0.35 

0.34 

± 

0.24 

2.03 

± 

2.28 

1.8 

38.7 ± 8.7t* 


STACK 

- 

3.66 

0.81 

0.85 ± 0.35 

- 

4.14 ± 0.71 

- 
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ZP-CDPS-261 

53.0826530 

-27.8664989 

3.40 

23.2 

24.2 

24.5 

± 

0.1 

27.1 

1.07 

0.40 

0.61 

± 

0.14 

3.54 

± 

0.80 

0.62 

± 

0.06 

1.21 

± 

0.25 

1.9 

12.1 ± 4.4t 

ZP-CDPS-400 

53.1025696 

-27.8606110 

4.10 

24.3 

25.1 

25.1 

± 

0.2 

23.9 

0.52 

0.33 

0.24 

± 

0.13 

1.45 

± 

0.78 

0.78 

± 

0.11 

3.40 

± 

1.40 

0.9 

31.3 ± 3.6t* 

ZP-CDPS-509 

53.1167717 

-27.8559704 

3.95 

24.2 

25.1 

25.0 

± 

0.0 

29.1 

0.41 

0.25 

0.31 

± 

0.06 

1.55 

± 

0.32 

0.52 

± 

0.05 

0.51 

± 

0.17 

1.0 

-4.5 ± 4.1 

ZP-COSMOS-12141 

150.0815277 

2.3637166 

4.00 

24.0 

24.7 

24.1 

± 

0.2 

18.8 

0.45 

0.34 

0.81 

± 

0.27 

3.58 

± 

1.09 

0.40 

± 

0.10 

4.92 

± 

1.35 

1.1 

0.9 ± 8.0 

ZP-COSMOS-3784 

150.1817627 

2.2390490 

3.58 

22.9 

23.9 

23.8 


0.1 

26.6 

0.36 

0.38 

0.53 


0.13 

3.40 

± 

0.78 

0.77 

± 

0.10 

1.88 

± 

0.33 

0.5 

-2.4 ± 10.2 

ZP-UDS-11279 

34.3843269 

-5.1402941 

3.72 

25.0 

26.6 

26.4 

± 

0.3 

4.5 

0.46 

0.32 

0.15 

± 

0.10 

0.96 

± 

0.54 

0.81 

± 

0.23 

1.00 



2.2 

29.3 ± 12.5 

ZP-UDS-4432 

34.3581772 

-5.2409291 

3.76 

23.8 

24.5 

24.2 

± 

0.2 

17.5 

0.83 

0.37 

0.75 

± 

0.39 

3.61 

± 

1.74 

0.46 

± 

0.11 

4.27 

± 

1.65 

1.5 

669.0 ± 10.7* 

ZP-UDS-4449 

34.3409157 

-5.2405076 

3.84 

23.1 

24.4 

24.9 

± 

0.1 

17.2 

0.41 

0.35 

0.44 

± 

0.10 

1.90 

± 

0.41 

0.38 

± 

0.07 

0.23 

± 

0.14 

1.0 

- 

ZP-UDS-4462 

34.3408661 

-5.2402906 

3.92 

23.0 

24.0 

24.0 

± 

0.1 

27.9 

0.39 

0.26 

0.39 

± 

0.09 

2.09 

± 

0.45 

0.60 

± 

0.08 

1.69 

± 

0.27 

0.8 

22.6 ± 9.4 

ZP-UDS-5617 

34.3407745 

-5.2240300 

4.17 

24.5 

26.0 

24.5 

± 

0.3 

6.3 

0.42 

0.37 

2.33 

± 

0.72 

10.74 ± 3.30 

0.45 

± 

0.18 

4.92 

± 

1.51 

1.3 

9.5 ± 9.7 

ZP-UDS-8379 

34.4104004 

-5.1821156 

3.77 

23.8 

25.2 

25.2 

± 

0.1 

14.0 

0.65 

0.25 

0.30 

± 

0.07 

1.50 

± 

0.34 

0.50 

± 

0.09 

0.52 

± 

0.28 

2.6 

355.8 ± 25.0* 

ZP-UDS-8399 

34.4105759 

-5.1825032 

3.44 

24.4 

25.3 

25.0 

± 

0.1 

11.9 

0.43 

0.23 

0.69 

± 

0.16 

2.28 

± 

0.49 

0.20 

± 

0.05 

0.14 

± 

0.17 

2.5 

106.6 ± 25.1* 

ZP-UDS-8580 

34.3544159 

-5.1797152 

4.07 

23.7 

24.6 

24.7 

± 

0.1 

19.8 

0.66 

0.26 

0.36 

± 

0.08 

1.82 

± 

0.37 

0.54 

± 

0.05 

0.18 

± 

0.09 

1.1 

7.1 ± 8.4 

ZP-UDS-9165 

34.3225441 

-5.1713767 

4.06 

23.4 

24.2 

24.6 


0.1 

33.8 

0.68 

0.31 

0.11 

± 

0.03 

0.66 

± 

0.14 

0.72 


0.09 

1.00 



0.3 

43.3 ± 10.1* 

STACK 

- 

- 

3.84 

- 

- 

- 
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Properties of 16 quiescent and 14 star forming galaxies 
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Fig. 1.— Left: Example galaxies at 2 4 of varying magnitude. Second: stacks of the quiescent and star forming subsamples, with the 

corresponding best-fit models and residuals after subtracting the models. Third: F814W/F125W/F160W stack color composites. Right: 
Stellar mass surface density profiles. Thin orange and blue lines represent individual measurements of quiescent and star-forming galaxies, 
respectively. Thick lines represent the stacks. The inset shows the surface brightness profiles of the stacks, with horizontal lines indicating 
So- brightness limits of 28.3 mag/arcsec^, measured in annuli of 0.06'' (0.43kpc) width at r > 1". The background limit for individual 
galaxies (dotted line) is 26.8 mag/arcsec^. 
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sources were masked. The final stacks were obtained by 
calculating the mean value at each pixel location of the 
image stamps. 

We ran GALFIT using the same input parameters as 
for the individual galaxies. Errors were estimated by 
bootstrapping, i.e. randomly selecting galaxies, recreat¬ 
ing the image stacks and rerunning GALFIT. 

In Figure [T] we show the stacks and examples of indi¬ 
vidual galaxies. The stack of quiescent galaxies is red¬ 
der than the stack of star-forming galaxies and has a 
more compact morphology. We also show stellar mass 
surface density profiles (E(M 0 /kpc^) = M{< r)/( 7 rr^)), 
obtained from the light profile measured in concentric 
apertures of radius r and assnming a constant mass-to- 
light ratio. For the stacked profiles we used the mean 
mass of the galaxies in each stack. They are consistent 
with the individual profiles within the uncertainties, sug¬ 
gesting that the stack does not reveal an extended low 
surface brightness component, down to a surface bright¬ 
ness limit of 28.3 mag/arcsec^. 

3.3. Contamination by AGN 

A substantial fraction of sources show signs of AGN 
activity either from X—ray detections or strong 24/rm 
(rest-frame 5/im, tracing hot dust). As WFG3/F160W 
(A = 1.5396/rm) corresponds to rest-frame 2960 — 3500A 
(UV) at 3.4 < z < 4.2, it could be that an AGN is 
dominating their central light, leading to small sizes of 
the single Sersic fits. 

In the quiescent sample we find four X—ray detected 
galaxies, two of which are spectr oscopically confirmed 
type-II QSOs (|Szokolv et al.ll200^ . Another has strong 
24/im, which could either point towards dust-obscured 
star-formation or AGN activity. Several have small pos¬ 
itive residuals after subtracting the best fit, suggesting 
the presence of a central point source. These 5/16 (31%) 
galaxies were re-fit with two components, a Sersic model 
and a pointsource-like model (represented by a Gaussian 
with FWHM= O.lpixels) to trace possible AGN light. In 
these models, the point source accounts for 4.3 — 68 % 
of the total light (with 57% and 68 % for the type-II 
QSOs, but on average 6 . 2 % for the remaining 3 AGN 
candidates). The average size of the Sersic component 
increases by 1.5 x (from a median Vf. = 0.13 ± 0.12" to 
re = 0.20 ±0.03"). 

Amongst the star-forming galaxies two are X—ray de¬ 
tected, and four are very bright at 24/im (L> 7 x 
IO^^Lq or SFR > 12 OOM 0 /yr). Re-fitting with a two- 
component model attributes 0.9 — 39.4% of the light to 
a point source, while the extended component changes 
in size by 0.65x (from Ve = 0.31 ± 0.15" to rg = 0.19 ± 
0.02"). We note that for the most extended sources, 
adding central light reduces the Sersic index ngersic of 
the extended component, and can result in a smaller rg. 

We additionally estimated the possible AGN contribu¬ 
tion from the galaxy SEDs. We first determine the best 
fitting powerla w bluewards of res t-frame 0.35/im and at 
observed 8/im (|Kriek et al.l 120091) . Then we fit the sum 
of the powerlaw and the original best-fit EAZY template 
to the data. The contribution of the AGN powerlaw tem¬ 
plate to F160W is 1.1 — 7.4% for the 5 quiescent galaxies 
and 0.9 — 2.9% for the 6 star-forming galaxies. 

While the two-component fits and SEDs indicate that 


a point source contribution is probably small, the true 
contribution and its effect on the sizes remain unclear. 

4. RESULTS 

We show the effective radius as a function of stellar 
mass in Figure [2J Quiescent galaxies at z ^ 4 are very 
compact, with a bootstrapped median size fg = 0.63 ± 
0.18kpc. When we remove AGN we find a similar result: 
rg = 0.57±0.18kpc. 

Star-forming galaxies have rg = 2.0 ± 0.60kpc. They 
are 3.2 ± 1.3 x larger than quiescent galaxies. Both sam¬ 
ples have a large spread in size, with some almost as large 
as at z ~ 0, showing that at z 4 the population is al¬ 
ready very div erse. On average the sizes lie well below the 
z ^ 0 relation (iMosleh et al.ll20l3l . by 6.0±1.7x for qui¬ 
escent and 2.4±0.7xfor star-forming galaxies. Quiescent 
galaxies are also 1.9 ± 0.7xsmaller than at 2 < z < 2.2. 


In Figure [3] we show Sersic index versus size for the 
z ~ 4 galaxies, and a sample at similar mass at 2 < 
z < 2.2. Star-forming galaxies have smaller Sersic in¬ 
dex, with on average Usersic = 1-3 ± 0.7, compared to 
naersic = 3.2 ± 1.2 for qifiescent galaxies. The difference 
between the two populations is also clear from the stellar 
mass density profiles in Figure [U with quiescent galaxies 
having steeper profiles and more centralized flux. In Fig¬ 
ure [3]we also plot (jAjmax, defined as the average stellar 
mass density inside the radius where S( Mf 7 ^/kpc^) falls of 
by a factor of two (|HoDkins et al.l[2010ll . with uncertain¬ 
ties from the Monte Garlo procedure described in section 

o 

Quiescent galaxies at z ^ 4 have a median {'S)max = 
3.3 ± 1.1 X lO^°M0kpc“^, much higher (^ lOx) than 
for star-forming galaxies: {^)max = 0.3 ± 0.1 x 
lO^°M 0 kpc“^, and more similar to 2 < z < 2.2 quiescent 
galaxies: (T,)max = 1-7 ± 0.3 x lO^°M 0 kpc“^. 

When stacking we find rg = 0.85 ± 0.35kpc (quies¬ 
cent) and Tg = 2.6 ± 1.2kpc (star-forming), and Sersic 
indices Usersic = 4.17 ± 0.90 and Usersic = 2.18 ± 2.03, 
respectively. The effective radius of the quiescent stack 
is slightly larger than the median of the individual galax¬ 
ies, by 1.3 ± 0.3X at < Icr significance, but overall the 
results are consistent. 

In Figure m we show the median sizes at the respective 
mean redshifts of the two subsamples. Gomparing with 
lower redshift, they continue to follow a trend of decreas¬ 
ing size with increasing redshift. Our control sample of 
galaxies at 2 < z < 3.4 with 10.5 < /McA < 11 

corres ponds well with the resnlts of Ivan der Wei et ahl 
(|2014f) . which suggest the same trend. 

We fit a relation of the form rg = A(1 ± z)^kpc at 
0 < z < 4, using the measurements of ivan der Wei et al.l 
(IMl at z < 2. We find rg = 5.08 ± 0.28(1 ± 
z)“i-44±o.08]jp(, £pj. quiescent and rg = 6.02 ± 0.28(1 ± 
2^-o.72±o.05]jpp star-forming galaxies. We note 
that our sample at z ^ 4 includes higher mass 

(logio(M/M 0 ) > 11) galaxies. If we remove the most 
massive galaxies, we find the same evolutionary relation. 

To test for incompleteness for diffuse galaxies, we red- 


















6 


QUIESCENT STAR-FORMING 



10.2 10.4 10.6 10.8 11.0 11.2 11 . 410.2 10.4 10.6 10.8 11.0 11.2 11.4 
log[M/Mg] log[M/Mg] 


Fig. 2. — Circularized effective radiu s for galaxies at 2: ^ 4. In purple and green we show our control sample at 2 < 2: < 2.6 and 
2.6 < 2 : < 3 and in orange the median of ivan der Wei et al.l (120141) at 2 < 2 : < 2.5. The black solid line is the z 0 relation of IMosleh et al.l 
1 ( 2013 ). X—ray detections and bright 24/im sources are indicated with stars and open circles. The median sizes are re = 0.63 zb O.lSkpc 
(quiescent galaxies) and Ve = 2.0 di 0.60kpc (star-forming galaxies). 

shift a stellar-mass matched sample with r > 2kpc and 
risersic < 2.5 at z ~ 2.5 to z = 3.7 and find 70% com¬ 
pleteness. 

5. DISCUSSION 

Our results show that the galaxies at z ^ 4 in 
this study obey similar relations between size and star¬ 
forming activity as galaxies at lower redshift. Quies¬ 
cent galaxies are compact, while star-forming galaxies are 
more extended and diffuse. The difference is also clear 
when selecting purely on size: if we define compactnes s 
as < 1.5 (|van der Wei et alll2014D . 

13/14 (93%) of massive compact galaxies would be clas¬ 
sified as quiescent, and 13/16 (81%) of larger galaxies as 
star-forming (Figure [3]). 

The number density of compact, /Mq) > 

10.55, quiescent galaxies at z 4 is 1.8 ± 0.8 x 
10“®Mpc“^, increasing by > 5x between 3.4 < z < 4.2 
and 2 < z < 2.2, towards 1.0 ± 0.3 x 10“"^Mpc“^. This 
suggests we are probing a key era of their formation, and 
we would expect to see their star-forming progenitors in 
abundance. 

Small effective radii for star-forming gal a xies have been 
repor ted at z = 2 — 3 (jBarro et al.ll2014J bl: (Nelson et al.l 
(20141) . They are rare in our sample: we find 1/14 with 
re/(M/lO^^M0)°-^^ < 1.5. On average, star-forming 
galaxies at z ~ 4 are twice as large as quiescent galaxies 
at z 2. If they are the direct progenitors of z < 4 
compact quiescent galaxies, we expect them to be simi¬ 
lar, not only in s ize, but also in Ser sic index and central 
surface density ([Nelson et al.((20l3 ). However, we find 
smaller nsersic for star-forming galaxies, while the cen¬ 
tral densities indicate that they must increase in (^)max 
by 5 — 10 X, to match the more cuspy profiles of z = 2 — 4 
quiescent galaxies. 


In a recent simulation, (Wellons et al.l (Ilustris; (20141) 
trace the evolution of galaxies to z = 2. They indeed 
identified two theoretical formation tracks: one in which 
a brief and intense central starburst prompted by a gas- 
rich major merger causes the galaxies’ half-mass radius 
to decrease dramatically. The second is that of a more 
gradual but early formation, with small galaxy sizes due 
to the higher density of the universe. In the second case, 
nearly all of the stellar mass is in place at z > 4. 

Comparing with the observations, we find that 19/44 of 
massive z ~ 4 galaxies are classified as quiescent, whereas 
all similarly massive galaxies in Illustris are still actively 
star-forming, with a typical SFR = 100 — 2OOM0/yr. 
This level of star-formation is ruled out at > 2>a by 
Herschel observations of the z ^ 4 quiescent galaxies 
(jStraatman et al.l(20141 ). At the same time, the fraction 
of compact galaxies in our sample is 47% ,versus 20% 
in Illustris. Hence massive galaxies appear to quench 
their star-formation earlier and to be more compact than 
in simulations. 

The paucity of compact star-forming galaxies at z ^ 4 
and their large median rest-frame UV size is puzzling. 
At face value it suggests that the rapid increase in num¬ 
ber density of compact quiescent galaxies cannot be ex¬ 
plained by simple shutdown of star-formation in typical 
star-forming galaxies of similar stellar mass. A possible 
solution is a rapidly forming dense core, i.e. a central 
starburst. Then the chance to observe the progenitors 
in our sample is small, as it is proportional to the du¬ 
ration of the main star-forming episode. For example, if 
compact cores of 2 < z < 2.2 quiescent galaxies formed 
at random times between 2.5 < z < 6, with a typical 
lOOMyr central starburst duration, their predicted num¬ 
ber density at z ^ 4 would be ^ 6 x 10“®Mpc“^. The 
observed number density of compact star-forming galax- 
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Stellar mass versus size. Bottom-left: Sersic index versus size. B ottom-right: stellar mass versus maximum stellar mass density. The 
horizontal dashed line is the empirical limit of lHookins et all 1)20101) . Only one z ^ 4 star-forming galaxy is compact. On average quiescent 
galaxies have smaller sizes, higher Sersic indices and higher central densities than star-forming galaxies. 


ies is 1.4 ± 1.4 X 10“®Mpc“^: smaller, but in a similar 
range given the large uncertainties. 

We note that the remarkably high fraction of quies¬ 
cent galaxies at z ^ 4 (Figure S]) is still uncertain. Cur¬ 
rent limits on the average dust-obscured SFR are weak 
(< l^Mr^lvrda). IStraatman et al.l[20l^ . hence some 
of the quiescent galaxies could be star-forming. Cosmic 
variance is significant 30%). Highly obscured massive 
star-forming galaxies might also be missed by nea r-IR 
surveys (e.g. IPaddi et al. 1 120091 : iCanuti et al.lfeoi^ . al¬ 
though the abundance and redshift distribution of such 
galaxies is still very uncertain. Finally, extended {r > 
3kpc) galaxies with small Usersic and low surface bright¬ 
ness are more difficult to detect than compact galaxies 
(e.g. Trujillo et al. 2006). 

We caution that the light prohles measured here may 
not be representative of the stellar mass distribution due 
to color gradients, with rest-frame UV sizes larger than 
rest-frame optical sizes. This would imply that the size 


evolution is stronger. However, using control sample at 
z ^ 3, we find no difference between U V and optical, 
consistent with Ivan der Wei et ahl l|2014f) , who show this 
effect is < 10% at z ~ 2 and decreasing with redshift. 

Galaxy sizes may also be overestimated if dust is ob¬ 
scuring a central starburst. Submm sizes of obscured 
starbursting galaxies could be small: < Ikpc (e.g. 
Ilkarashi et al.ll20i^ ISimpson et al.l[20T5[l . A direct com¬ 
parison of ALMA submm and rest-frame optical/UV 
morphologies for the same objects with measured stel¬ 
lar mass will reveal the effect of dust obscuration on 
UV/optically measured galaxy sizes. 
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